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Abstract. The case for substantial far infrared ice emission in local ultra- 
luminous infrared galaxies, expected based on the presence of mid-infrared ice 
absorption in their spectra and the known far infrared optical properties of ice, is 
still largely unsupported by direct observation owing to insufficient far infrared 
^ | spectral coverage. Some marginal supportive evidence is presented here. 

, A clear consequence of far infrared ice emission is the need to extend the 

l/^ ■ range of redshifts considered for submillimeter sources. This is demonstrated 

via the example of HDF 850.1. 

The solid phase of the ISM during reionization may be dominated by ice, 
and this could lead to the presence of reionization sources in submillimeter source 
catalogs. Submillimeter sources not detected at 24 fim in the GOODS-N field are 
examined. Two candidate reionization sources are identified at 3.6 //m through 
possible Gunn-Peterson saturation in the Z band. 
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. 1. Ice emission 
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. The discovery of ice absorption features in the spectra of ultraluminous infrared 

galaxies (ULIRGs; Spoon et al. 2002; Imanishi & Maloney 2003; Imanishi, 
Dudley & Maloney 2006) implies that ice is emitting in the far infrared (FIR). 
Qh| The obervational consequences are determined by the portion of FIR emission 

q ■ that is owing to mid-infrared absorption (often the majority) and the ratio of 

mantle-to-core volumes of the emitting dust as well as the FIR optical depth. 
For ULIRG FIR dust temperatures one expects the broad 150 fim phonon mode 
(Berti et al. 1969; Curtis et al. 2005) of ice to modify the FIR spectral energy 
distribution (SED) through emission. An example is given in Fig. la. 
\ This changes at the higher dust temperatures required at higher redshift 

where the CMB temperature exceeds (optically thin) ULIRG temperatures. 
Then, shorter wavelength ice features are more important to SED modification. 
Fig. lb. shows HDF 850.1, the brightest submillimeter source in the Hubble Deep 
Field, plotted against an ice emission model. The redshift, z — 12.6, is 13 con- 
fidence intervals beyond the estimate of Dunlop et al. (2004): z = 4.1. Thus, 
ice emission could greatly broaden the range of redshifts that need to be consid- 
ered for submillimeter sources. Fig. 3a. shows the expected 850 /im brightness 
of sources with the luminosity of ULIRG IRAS 00188-0856 if they are similar 
to Arp 220 (long dashed line), have substantial ice and a temperature of 30 K 
(solid line), or are similar to the model shown in Fig. lb. (50 K; short dashed 
line). It is notable that the last model selects for z = 13. 

2. What is the early ISM like? 

When metals are first released into the interstellar medium the oxygen-to-carbon 
ratio is high and the 56 Ni production is frequently low. For a 186 M progenitor 
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Fig. la. Two fits for the FIR SED of ULIRG 
IRAS 14348-1447 are shown. The dashed line 
reproduces the fit given by Klaas et al. 2001 
and has a 100 /im optical depth of 5 and (3 — 2 
while the solid line is an ice mantle-core emis- 
sion model taken from Aannestad (1975). The 
broad feature centered near 150 /xm is hinted 
at in the ISOPHOT data though one fit is not 
better than the other. This source shows 6 /im 
absorption e.g. Charmandaris et al. (2002). 
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Fig. lb. Data from Dunlop et al. (2004) plus 
a new Z band limit (Dudley et al. 2006; filled 
circles; lensing corrected) are plotted together 
with new Spitzer GOODS data for the lensing 
galaxy (filled diamonds; relative scale). Solid 
lines show an ice emission model and a star- 
burst model (Leitherer et al. 1999). Dashed 
lines and the shaded region represent evolu- 
tion of the lensing eliptical galaxy e.g. Athey 
et al. (2002). 



(90 M core) the yield ratios are O:Si:C:S: 56 Ni=123:29:15:10:l by number (Heger 
& Woosley 2002) so that the onset of dust formation may be rapid and occur in 
an oxygen rich environment. If grain temperature is determined by kinetic and 
chemical interactions when the mainly oxygen eject a mix with the surrounding 
mainly hydrogen atmosphere (Baraffe, Heger & Woosley 2001) for example, 
rather than by ongoing radioactivity, grain surface and gas phase chemistry 
might lead to the condensation of ice mantles at temperatures high enough to 
lead to the formation of crystalline rather than amorphous ice such as is observed 
to occur in some evolved stars. 

If all or most oxygen enters the solid state, the volume of ice at this early 
phase would be larger than the volume of more refractory material. Further, 
if the injection of such material is needed to facilitate the formation of the 
generation of stars that is proposed to be 4 Gyr old at z =1.5 (Jimenez et al. 
2000) or >2 Gyr old at z =2.5 (Stockton et al. 2004) then the formation interval, 
if set by this initial enrichment, may be brief enough that a substantial fraction 
of stars are found simultaneously within their natal envelopes which may yet 
contain crystalline ice. 



3. Does SCUBA see to reionization? 

The positional error cicles shown in Fig. 2 for SCUBA 850 /im sources cataloged 
by Wang et al. (2004) are remarkable because so many of them are empty at 24 
/im. Six are definitely clear of obvious sources and an additional four (sources 27, 
28, 41 and 44) are borderline. If every SCUBA source has a 24 /im counterpart, 
we expect this situation to occur 2.25 times using 2a radii error circles. And, 
that is in the absence of coincidental association. So, some number greater than 
zero of the 850 /im sources may not have 24 /im counterparts. 

At the sensitivity of the GOODS MIPS image that implies z > 4 using 
ULIRG templates. Chapman et al. (2005) argue that radio non-detected sub- 
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Fig. 2. Spitzer GOODS data for 
the 45 Wang et al. (2004) HDF- 
GOODS submillimeter sources. 
The darker red-numbered bands 
are MIPS 24 /im images, and 
below these are the correspond- 
ing 8, 5.8, 4.5 and 3.6 /xm IRAC 
data. The numbers are as given 
by Wang et al. with their 
(green) 2a radii positional er- 
ror circles. Dots (pink) are ra- 
dio positions of claimed (Chap- 
man et al. 2005; Wang et al.) 
associations. HDFGOODS 850 
/xm sources 8, 14, 18, 21, 35 
and 45 do not have obvious 24 
/xm couterparts within the er- 
ror circles. Yellow arrows in- 
dicate sources we have discov- 
ered at 3.6 /xm that do not have 
Z band counterparts, similar to 
HDF 850.1 



millimeter sources have a similar redshift distribution to those radio detected 
sources for which they have obtained optical spectra. However, their argument 
depends on the SEDs of the sources being similar to normal galaxies, and Fig. 
3b. shows that this may not be the case. So, the MIPS data together with a 
reanalysis of sources with known redshifts suggest a higher redshift population 
could be present in submillimeter catalogs. 






Fig. 3a. 850 /im brightness is plotted against 
redshift for a typical SED model (long dashed 
line; Scoville et al. 1991), an ice emission 
model for IRAS 00188-0856 (solid line; Dud- 
ley et al. 2006) and the model shown in Fig. 
lb. The right hand scale accounts for lensing 
amplification for HDF 850.1. The arrows are 
redshift estimates for this source from Dunlop 
et al. and Dudley et al. 



Fig. 3b. 850 /xm-to-radio ratios are caclulated 
for each of the listed ULIRGs at the redshifts 
of the sources cataloged by Chapman et al. 
(2005) for 1 < z < 3.2. The Chapman et 
al. sources are assigned to bins based on the 
ULIRG distribution. The histogram is shown 
by the solid line. The dashed line shows a 
bin width correction. Bin average redshifts are 
listed above the lower axis. 

We find sources in the 24 /im blank fields for HDFGOODS 850-8 and 850-14 
which are similar to HDF 850.1 in that they are detected in the "NIR" (3.6 /im) 
but not at / or Z band (Fig. 4). I — [3.6] > 5 and 4.5 (Vega) respectively. 
Should ice play the role suggested in Fig. lb. in these sources, making them 
detectable at 850 fim but not too luminous to be precursors of the known old 
elipticals at z =1.5 and 2.5, then these could be reionization sources. If Gunn- 
Peterson saturation explains the Z band non-detections then the redshifts of 
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Fig. 4. Spitzer IRAC 3.6 fim images for the fields of HDFGOODS 850-8 and 850-14 are 
reproduced from Fig. 2 together with Hubble ACS Z and / band images. The 3.6 /urn sources 
indicated with the yellow arrows are absent at the shorter wavelengths. Green contours are 
3.6 /im data. 

these sources would be greater than 7 and thus larger that the z ~6 onset of the 
Gunn-Peterson trough (Fan et al. 2004) while the redshift of last star formation 
is z ~13 for the old ellipticals in a WMAP cosmology (Spergel et al. 2003). 
This coincides with the redshift range for reionization estimated by Kogut et al. 
(2003). These sources may prove easier to observe than HDF 850.1 owing to the 
absence of a lens. 
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